Coulomb screening and collective excitations in biased bilayer graphene 

Xue-Feng Wan£] 

Department of physics, Soochow University, 1 Shizi Street, Suzhou, China 215006 



o 

(N 

a 
in 



c 

o 
o 



> 

(N 

o 
o 



X 



Tapash Chakraborty 

Department of Physics and Astronomy, The University of Manitoba, Winnipeg, Canada, R3T 2N2 

We have investigated the Coulomb screening properties and plasmon spectrum in a bilayer 
graphene under a perpendicular electric bias. The bias voltage applied between the two graphene 
layers opens a gap in the single particle energy spectrum and modifies the many-body correlations 
and collective excitations. The energy gap can soften the plasmon modes and lead to a crossover 
of the plasmons from a Landau damped mode to being undamped. Plasmon modes of long lifetime 
may be observable in experiments and may have potentials for device applications. 
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Bilayer graphene (BLG) has attracted much attention 
due to its unique electronic characteristics, distinct from 
the Dirac gas in monolayer graphene and the Fermi gas in 
traditional semiconductor quantum wells [H-Q- I n addi- 
tion, an energy gap between the conduction and valence 
bands of a BLG can be opened and tuned by introducing 
an electrostatic potential bias between the two graphene 
layers 0-[ll|. This can be easily realized via one or more 
external gates to perpendicularly bias BLG and make 
it a potential component for integrated electronics. It 
is then very intriguing to understand some fundamental 
properties such as correlation and screening properties of 
electron gases in a biased BLG. As collective excitations, 
plasmon modes are a direct result of electronic correlation 
due to Coulomb interaction between electrons. Experi- 
mental detection of plasmon modes has recently become 
feasible and has been used to determine the dynamical 
behavior of electrons in graphene layers (l2l - ll4j |. 



Previously, assuming zero or non-zero spin-orbit inter- 
action induced energy gap, we have studied the Coulomb 
screening and collective excitation spectrum of intrin- 
sic and doped monolayer graphenes at zero and finite 
temperatures in the random phase approximation (RPA) 
[l5l | . Later, Qaiumzadeh and Asgari [l6| assumed an un- 
specified energy gap of arbitrary width for doped mono- 
layer graphene and studied the corresponding ground- 
state properties at zero temperature in RPA. They con- 
cluded that the conductance and charge compressibility 
decrease with the band gap. Furthermore, a THz source 
has been proposed based on the stimulated plasmon emis- 
sion in graphene [l8| and the absorption of THz electro- 
magnetic radiation in gapped graphene has been esti- 
mated 19]. On the other hand, the Coulomb screening 
and the collective excitations in zero gap BLG have been 
studied in our previous work at zero and finite tempera- 
tures @ and by Hwang and Das Sarma Q for the zero 
temperature case. In this paper we report on our studies 



of the correlations, screening, and the plamson spectrum 
of electron gases in a biased BLG. 

In the effective-mass approximation the Hamilto- 
nian describing electrons of moderate energies in the K 
valley of a biased BLG reads 
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with k± = k x ±ik y = — iV^^zVy and k = (k x ,k y ) being 
measured from the K point. The effective mass of the 
quadratic term is m * = 2fi 2 7i/(3a 7o) 2 ~ 0.033m with 
too the free electron mass, clq — 1.42 A the C-C bond 
length on the graphene layer, 70 = 3.16 eV the intra-layer 
coupling, and 71 = 0.4 eV the direct inter-layer coupling. 
The second term arises from the electrostatic potential 
bias U between the two graphene layers separated by a 
distance d = 3.35 A. The Hamiltonian is obtained by 
keeping only the linear term of the Tayler expansion on 
the small energy value in unit of 71, so it is valid for 
electrons of energy less than 0.4 eV which is adequate in 
our case. The indirect inter-layer coupling is neglected 
since it affects only the energy band in the range of less 
than 2 meV from the middle of the conduction-valence 
band gap Q. 

The eigenenergy of the above Hamiltonian is 

with the eigenfunc- 



tions "ffc 1 ^) 
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Here is the azimuth of the vector k, i.e., tan 6k = k y /k x , 
and ak indicates the ratio of the kinetic energy to the po- 
tential bias with tana*; = h 2 k 2 /(m*U). The conduction 
band which touches the valence band at k = in unbi- 
ased BLG becomes separated from it by an energy gap 
equal to the potential bias U. This gap converts the BLG 
from a semimetal into a semiconductor and accordingly 
modifies the optical and electric properties of the elec- 
trons inside. For finite U, the density of states of the 
BLG diverges on the edge of the energy gap = {7/2. 
At zero temperature, the carrier density N in a BLG of 
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Fermi energy E F is N = ±^f^/E 2 F - U 2 /4. 

Following the well-established formalism for spin sys- 
tems [20| . we obtain the dielectric matrix of a biased BLG 
in the form of a unit matrix multiplied by a dielectric 
function e(q,co) = 1 — v q Tl(q, uj) with the bare Coulomb 
interaction v q = e 2 /(2eoq) and the electron-hole propa- 
gator 



n(q J w)=4 



E 

A.A'.fc 



!^' v (<?)l 2 



(2) 



The factor four comes from the degenerate two spins 
and two valleys at K and K' , f(x) is the Fermi 

function, and the vertex factor reads \g^' X (<z)| 2 = 
i [1 + AA' cos «fc cos otk+q + XX' sin a k sin a^+q cos(2t9fc — 

29 k+q )}. At q = or q = -2k, \g£ x ' (q)\ 2 - (1 + AA')/2. 
Similar to unbiased BLG, the interband vertical and back 
scatterings are both forbidden but the intraband back 
scattering is allowed in biased BLG. 

It has been shown that the interlayer indirect C-C in- 
teraction introduces anisotropic fine structures near the 
Fermi energy in the range of 2 meV and leads to some 
interesting dielectric and collective phenomena 0. For 
systems with energy gap U > 5 meV or with Fermi en- 
ergy Ep satisfying \Ep — kT\ > 3 meV, this anisotropy 
becomes negligible. For large U comparable to vi, the 
effect of the "Maxican hat" at the bottom (top) of the 
conduction (valence) band [1] should be taken into ac- 
count. Nevertheless, for moderate U and Ep, the model 
described here should be valid. Furthermore, we assume 
that the BLG is far enough from the substrate and the 
gate so a unit background dielectric constant is used in 
the calculation. 

In intrinsic BLG where no net carrier exists, i.e., N = 
and the Fermi energy Ep — 0, intraband scattering is 
only allowed at non-zero temperatures. In Fig. [I] we have 
shown that the real part (e r , solid curve) and imaginary 
part (si, dotted) of the dielectric function versus the en- 
ergy in an intrinsic system with potential bias U = 5 
meV at a finite temperature 77 K for (a) q — 0.005 x 10 8 



m" 1 and (b) q = 0.5 x 10 8 



In the small q case, the 



intraband scattering introduces a dip for e r and a peak 
for £j of low energy as illustrated in the insets. Con- 
sequently, there exist two plasmon modes, one Landau 
damped and one almost undamped. The depth and the 
width of this real part dip increase with the temperature 
indicating the increase of intraband scattering strength 
and also the energy of the undamped plasmon mode. At 
uj = \JlJ 2 + q 4 /lQm*, the threshold of interband single- 
particle excitation continuum (SPEC), £j steps up and a 
sharp peak of e r is observed thanks to the flat bottom 
and top of the energy bands. This peak may introduce 
additional plasmon modes and is similar to the case in 
gapped monolayer graphene [lBj]. As q increases, the e r 
dip (e, peak) due to intraband scattering shifts quickly 
to the higher energy side while the e r peak (ei step) due 
to the interband scattering moves only slowly. As a re- 
sult, the well separated intra- and inter-band structures 
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FIG. 1: e r (solid) and Si (dotted) are plotted versus uj in 
biased intrinsic BLG at temperature T=77K for (a) q — 
0.005 x 10 8 m" 1 and (b) q = 0.5 x 10 8 m" 1 . The poten- 
tial bias is U=5 meV and the Fermi energy Ef = 0. The 
details at low frequency for small q is shown in the inset of 
(a). 



O 
II 

3 



I 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 




j E 20 

- interband 



q(10 8 m" 



1< 



FIG. 2: The intra- and inter-band contributions to the real 
part of the zero temperature propagator function —II in doped 
BLG with carrier density N = 10 12 m~ 2 versus the wavevec- 
tor q are illustrated for U — (solid), 20 (dashed), 40 (dot- 
ted), and 60 (dash-dotted) meV. The corresponding energy 
bands for different U is shown in the inset. iVo = 2m* /n is 
the density of states of unbiased BLG. 



at small q mix with each other and then separate again 
when q increase as shown in Fig.[Hb). 

The effect of a bias on the zero-temperature propaga- 
tor — II(<7, uj) in doped BLG Q with a fixed carrier den- 
sity = 10 12 m~ 2 is studied in Fig.[5J The interband 
contribution decreases with the bias potential U as the 
energy gap widens. The intraband contribution, on the 
contrary, increases with the bias since the energy disper- 
sion leads to an enhancement of the density of states near 
the Fermi energy. Here one of the characteristics in BLG 
against in monolayer graphene 0, [1] is the strong back 
scattering of electrons on the Fermi surface which results 
in an intraband peak at q = 2kp — 3.54 x I0 8 m _1 . If 




FIG. 3: The plasmon spectrum (solid curves for undamped 
or slightly damped mode and dashed for Landau damped 
mode) and single particle continuum spectrum (light shadow 
for interband SPEC and dark shadow for intraband SPEC) 
at zero temperature under potential bias (a) U = 30 and 
(b) U = 60 meV are illustrated for BLG with TV = 10 12 
m~ 2 . The plasmon dispersion in the long wavelength limit, 
uio — \J (e 2 /47reo) / (Nirq/m*) is also presented. 

the Fermi energy remains fixed as the bias increase, the 
carrier density decreases and the intra- (inter-) band con- 
tribution at large (small) q becomes less sensitive to the 
bias and decreases (increases) with the bias. 

The plasmon modes are obtained by solving the zeros 
of the real part of the dielectric function e r (q,uj) = 
and the corresponding imaginary part represents the 
damping rate of the plasmon modes. In Fig.[3j we plot 
the typical spectrum of plasmon modes (solid and dashed 
curves) in doped BLG at zero temperature under po- 
tential bias (a) U — 30 meV and (b) U = 60 meV. 
The uppe r light shadow is th e interband SPEC ed ged 
at cj = yJU 2 + k%/m* 2 /2 + y/U 2 + (fc F - q) 4 /m* 2 and 
the low er dark shadow is th e int raband SPEC edged at 
u = y/XP + (k F + qY/m* 2 ~ y/U 2 + k 4 F /m* 2 /2. One 
undamped mode is located in the SPEC gap due to fi- 
nite Fermi energy. In the long wavelength limit, its dis- 
persion is th e same as that of Fer mi 2D gas of two val- 
ley, wo = y/ (e 2 /47reo)/ (Nirq/m*). However, compared 
to the plasmon dispersion in BLG without bias which 
is just slightly modified from that of Fermi 2D gas, this 
dispersion is greatly softened for finite q as also shown 
in Fig. 2] Our numerical analysis shows that this is a re- 
sult of the deformation near the bottom and top of the 
energy bands. Note that the lowered plasmon group ve- 
locity may be helpful for making a stimulated plasmon 
oscillator [l8| . A Landau damped mode is located just 
below the intraband SPEC edge as usually happens in 
traditional 2D Fermi gas but is pushed to lower energy 
at larger q. The undamped mode can enter into the in- 
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FIG. 4: The energy u> versus the potential bias U of zero 
temperature plasmon modes (solid for undamped and dashed 
for Landau damped) at fixed Ef = 36 meV (left panels) or at 
fixed N = 10 12 m~ 2 (right panels) for q = 0.05 in (a) and (d), 
0.5 in (b) and (e), and 1 x 10 8 m -1 in (c) and (f), respectively. 



terband SPEC and becomes a slightly damped mode in 
some cases as shown in Fig.^a) under U — 30 meV or 
merges with the damped mode and disappears near the 
cross of intra- and inter-band SPEC edges as shown in 
Fig.EKb) under U = 60 meV. 

With the plasmon spectrum in mind, we now explore 
how U affects the energy and damping properties of the 
modes. In the left panels of Fig.|U we show uj versus U 
at several typical q when keeping the Fermi energy con- 
stant. As in Fig. [31 the light shadow indicates the inter- 
band SPEC and the dark shadow for the intraband one. 
At small q as illustrated in (a), there is one undamped 
plasmon mode with energy located inside the SPEC gap 
of which the width is about 2Ep and one damped mode 
of low frequency. When U reaches and passes 2Ep, the 
Fermi level drops below the conduction bottom and the 
two plasmon modes merge and disappear. At larger q, 
the intraband SPEC edge shifts up and the interband one 
shifts down for U < 2Ep and u) increases as shown in (b) 
and also in Fig.|3] Then the two SPECs will merge and 
the previous undamped plasmon mode enter the inter- 
band SPEC and become slightly damped. In this case, 
we may open the SPEC gap again by applying a stronger 
bias and transfer the slightly damped plasmon mode into 
a undamped as shown in (c) . The uj versus U curve forms 
a shoulder when it meets the interband SPEC reflecting 
the strong coupling between the single particle and col- 
lective excitations as also shown in other cases @, [TH, [2(| ■ 

If N remains constant as shown in the right panels, 
the Fermi vector is also constant but the Ep shifts up 
with U. This is clearly shown in (d) by the interband 
SPEC edge of small q which is located near oj = 2Ep. 
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The undamped plasmon mode continues to exist as U 
increases and its energy varies slowly. This is because 
Ep is always higher than the conduction bottom with a 
constant Fermi vector, uj decreases with U as the effec- 
tive mass near Ep increases. For a large q the plasmon 
mode located inside the interband SPEC and is slightly 
damped at small U, one can always make it undamped 
by increasing the bias and widening the gap between the 
intra- and inter-band SPECs as shown in (f). 

When an external gate voltage is applied to a BLG, the 
carrier density varies with the gate voltage as well as the 
energy gap. |7Hig| Although N and U can be dependent 
on each other in a nontrivial way, our result suggests that 
u> is proportional to \/~N in almost the same way in both 
doped and undoped BLG. This happens because u is 
mainly determined by N as illustrated in the right panels 
of Fig. [4] The variation of U of small amount affects ui 
only in a very limited scale. Nevertheless, as shown in 
Fig.|4j the higher U opens an wider energy gap in the 
SPEC and prolongs the lifetime of the plasmon modes. 
In other words, a gate voltage can vary the imaginary 
part of the dielectric constant at the plasmon energy and 
the effect may be observed in experiments. 



In summary, a potential bias can be applied between 
the two graphene layers of a bilayer graphene with the 
help of a gate voltage. We have studied the effect of the 
potential bias on electronic correlations, Coulomb screen- 
ing, and collective excitations at both zero and finite tem- 
perature. The potential bias opens a gap in the single 
particle energy spectrum and makes the semimetal bi- 
layer graphene a semiconductor. As a result the dielectric 
function for the Coulomb interaction and the propagator 
function are modified significantly. The potential bias 
also opens a gap in the single-particle excitation spec- 
trum and softens the collective excitation modes. This 
may result in undamped collective excitation modes that 
are observable in experiments. In the single gate config- 
uration, the doping and gate voltage can vary the poten- 
tial bias and the carrier density of the bilayer graphene 
and manipulate the energy and lifetime of the collective 
excitation modes inside. 
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